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Abstract 

We present a novel method of simulating wave effects in graphics 
using ray-based renderers with a new function: the Wave BSDF 
(Bidirectional Scattering Distribution Function). Reflections from 
neighboring surface patches represented by local BSDFs are mu- 
tually independent. However, in many surfaces with wavelength- 
scale microstructures, interference and diffraction requires a joint 
analysis of reflected wavefronts from neighboring patches. We 
demonstrate a simple method to compute the BSDF for the entire 
microstructure, which can be used independently for each patch. 
This allows us to use traditional ray-based rendering pipelines to 
synthesize wave effects of light and sound. We exploit the Wigner 
Distribution Function (WDF) to create transmissive, reflective, and 
emissive BSDFs for various diffraction phenomena in a physically 
accurate way. In contrast to previous methods for computing in- 
terference, we circumvent the need to explicitly keep track of the 
phase of the wave by using BSDFs that include positive as well as 
negative coefficients. We describe and compare the theory in re- 
lation to well understood concepts in rendering and demonstrate a 
straightforward implementation. In conjunction with standard ray- 
tracers, such as PBRT, we demonstrate wave effects for a range 
of scenarios such as multi-bounce diffraction materials, holograms 
and reflection of high frequency surfaces. 

1 Introduction 

Diffraction is a common phenomenon in nature when dealing with 
small scale occluders. It can be observed on animals, such as feath- 
ers and butterfly wings, and man-made objects like rainbow holo- 
grams. In acoustics, the effect of diffraction is even more significant 
due to the much longer wavelength of sound. In order to simu- 
late effects such as interference and diffraction within a ray based 
framework, the phase of light or sound waves needs to be integrated 
into those methods. 

We introduce a novel method for creating Bidirectional Scattering 
Distribution Functions (BSDFs), which efficiently simulate diffrac- 
tion and interference in ray-based frameworks. The reflected or 
scattered radiance of a ray indirectly and independently encodes the 
mutual phase information among the rays, conveniently allowing 
for interference after multiple bounces. Our BSDFs, derived from 
the Wigner Distribution Function (WDF) in wave optics, abstract 
away the complexity of phase calculations. Traditional ray-based 
renderers, without modifications, can directly use these WBSDFs. 
The implementation of the WBSDFs does not require a reader to 
fully understand the theory of the WBSDF derivation. A reader 
who is interested solely in i mplementing our new shaders can pro- 
ceed directly to section [TT] Specific technical contributions are as 
follows: 

1 . A method to compute WBSDF from microstructures 

2. Formulation of the rendering equation to simulate diffraction 
and interference of light 

3. A practical way to support interference and diffraction in 
global illumination with importance sampling 

4. Compatibility demonstration with ray-b ased render ers by cre- 
ating a material plugin for PBRT | |Pharr and Humphreys 2004 J 



Figure 1: We generalize the rendering equation and the BSDF to 
simulate wave phenomena using an unmodified PBRT rendering 
framework. The new Wave BSDF behaves like a local scattering 
function, creates interference globally, and allows easy integration 
into traditional ray based methods. 



5. Application of the new model to simulate rainbow holograms 
and camera optics 

6. Comparison to other diffraction techniques and suggestions 
for further improvement 

7. An extension towards sound rendering (see appendix) 

1.1 Related Work 

Light Propagation in Optics: In wave optics, light is described as 
an electromagnetic field with amplitude and phase. The Huygens- 
Fresnel principle is often used to represe nt wave propaga tion, 
which is a convolution of point scatterers [Goodman 2005| . In 
contrast, geometrical optics treats light as a collection of rays. 
Among the extensive efforts to connect wave and ray optics | Wolf| 
1978], notable ones are the generalized radiance proposed by 
Walther 1 1973] and the Wigner Distribution Function |B astiaans| 
|1977| , where light is described in terms of local spatial frequency. 
Although the generalized radiance or the WDF can be negative, 
it exhibits convenient properties that explains diffraction rigor- 
ously | |Bastiaans 1997j . 

Traditional Light Propagation in Graphics: Ray-based ren- 
dering systems, e.g., ray tracing [Whitted 1980 |, are popular for 
rendering photorealistic images in computer graphics due to their 
simplicity and efficiency. They are particularly convenient for sim- 
ulating reflection and refraction. In addition, in combination with 
global illumination techniques such as photon mapping | Jensen] 
|1996|pensen and C hristensen 1998], caustics and indirect light can 
be constructed. The idea of using negative light in the rendering 
equation ha s been proposed for visibility calculations | |Dachsbacher| 
|et al. 2007) , but not for interference calculations. 

Wave-based Image Rendering: Moravec proposed using a 
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Table 1: Overview of different rendering strategies for diffraction and interference. Our new technique addresses all of the required features 
with high efficiency and simplicity. OPD is the only technique that encompasses all of the effects that WBSDF simulates, but WBSDF is more 
efficient and easier to use with standard rendering pipelines. 



wave model to render complex light transport efficiently [19811. 
Stam implemented a diffraction shader based on the Kirchhoff 
integral 1 1999 1 for random or periodic patterns. Other variations of 
diffraction based BRDFs were created for rendering specific types 
of materials |Sun et al. 2 000, Sun 2006 1. Some other examples 
are b ased on the Huygens-Fresnel principle [Lindsay and Agu| 
|2006| . These, however, all compute diffraction and interference 
for an incoming and outgoing direction instantaneously at the 
location of reflection, whereas we defer the calculations of 
interference to a later stage. Zhang and Levoy were the first to 
introduce the connection between rays and the WDF in computer 
graphics [2009 1. The Augmented light field (ALF), inspired by 
the WDF, was presented by Oh et al. [2010, 2009], describing 
how transmission of light through a mask can be modeled using 
ray based rendering techniques. Optical Path Differencing (OPD) 
techniques keep track of the distance a ray travels and calculate its 
phase. Ziegler et al. developed a wave-based framework |^2008J, 
where complex values can be assigned for occluders to account 



for phase effects. They also implemented hologram rendering 
based on wave propagation (with the spatial frequency) [2007 1. 
Edge diffraction allows speedup by first searching for diffracting 
edges a nd then creating new source s at those positions [Freniere 
let al. 1 999, Tsingos Mar ch 2000| . In contrast, our WBSDF 
indirectly encodes the phase information in the reflected radi- 
ances and directions by introducing negative real coefficients. 
Therefore, no modification to the rendering framework is necessary. 

Wave-based Audio Rendering: For efficient rendering of sound, 
many techniques assume that high frequency sound waves can 
be modeled as rays. The two main diffraction models used in 
these geometric simulations are the Uniform Theory of Diffraction 
(UTD) [Kouyoumjian and Pathak 1974, Tsingos et al. 2001 1 and 
the Biot-Tols toy-Medwin (BTM) model [Hoth ersall et al. 1991| 
[Torres et al. 2001 1. These techniques, which are often referred 
to as edge-diffraction, first search for diffracting edges and then 
create new sources at those positions. UTD based-methods scatter 



Figure 2: Relationship between spatial frequency and incident an- 
gle, (a) Spatial frequency u of incoming light is dependent on the 
incident angle Oi and the wavelength A of the light, (b) The steeper 
incoming angle, the higher spatial frequency becomes. 
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incoming sound beams in a cone around the edge, where the 
half-angle of the cone is determine d by the a ngle that the ray hits 
the edge [Chandak et al. 2008, R ick and M athar 2007 1. They 
are often preferred over BTM because of their efficiency, but are 
less accurate in lower frequencies. BTM is more accurate and i s 
also applicable for non-geometric acoustics [Torres et al. 2001) , 
but is not applicable for interactive renderings due to its complexity. 



An overview and comparison of several diffraction and interference 
simulation methods is presented in Table [T] and Section |5] 

2 Formulation of a Wave Based-BSDF 

Consider the rendering equation [Kajiya 1986| : 



Lo(x,^o, A) = Le(x,^o, A) + / p{x 



(1) 

where Lo (x^Oo, A) is the total outgoing light from the point x in the 
direction Oo with wavelength A. Le describes the light emitted at 
the point x in the same direction, and Li denotes the incoming light 
at the point from a certain direction Oi . The function p represents 
the proportion of scattered light for a given incoming and outgoing 
direction, and is often referred to as the Bidirectional Scattering 
Distribution Function (BSDF). 

As Eq. [T] only takes the intensity of rays into account but not 
phase information, interference between multiple rays cannot be 
described unless travel distances of individual rays are tracked. Pre- 
vious methods have described how to add local diffraction and in- 
terference effects to the Bidirectional Reflection Distribution Func- 
tion (BRDF) IS tam 1999 J . However, a traveling ray does not carry 
phase information and therefore is unable to interfere in a later 
stage. It is challenging to represent this deferred interference in ray 
space. The Wigner Distribution Function, however, is a convenient 
method to understand and represent diffraction and interference. 

The Wigner Distribution relates the spatial x and spatial frequency 
u content of a given function. While it can be applied to many dif- 
ferent functions (i.e. music, images, etc.), it is a particularly useful 
description of a wave. The Wigner Distribution of a ID complex 
function of space t{x) can be defined as 
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where * is the conjugate operator. The function 



dx^ (2) 
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Figure 3: Creating WBSDF from micro structure, (a) A simple si- 
nusoidal grating is represented as a phase function, (b) Its WDF in 
position-spatial frequency {x,u) domain, where red defines posi- 
tive values and blue negative (c) We convert the vertical slices into 
plots for incoming and outgoing angle, creating the WBSDF ( shown 
at three different locations on the microstructure) (d) Rendering us- 
ing this WBSDF. Top-right shows a photo for visual validation. 



is often called the mutual intensity j Bastiaans 2009| and is the 
correlation function of a complex-valued microstructure geometry 
t{x). Note that after the Fourier transform of the mutual intensity, 
the WDF contains only real values, positive as well as negative, 
since the mutual intensity is Hermitian. In this section we use a ID 
function as input to explain the concept of using the WDF, but the 
extension to a 2D input signal is straightforward as 



Wt 



where 



J{x,y,x,y)=t(x+—,y+ — 



(4) 



(5) 



The essence of this representation is that a complex wavefront is 
decomposed into a series of plane waves with spatial frequency u 
and a real-valued amplitude. This spatial freq uency corr esponds 
to the directionality of the parallel wavefront [Goodman 2005] as 
illustrated in Figure |2|a-b). The local spatial frequency is related to 
the wavelength and the direction of wave; which is normal to the 
wavefront, as 

sin Oi 

«=— . (6) 

The Wigner Distribution is often used for describing the complex- 
valued wavefront arising from surfaces like a grating. If a plane 
wave (i.e. light from a point light source at infinity) hits a sur- 
face, the outgoing wavefront function Rt can be described with this 
Wigner Distribution Function using Eq.[2]as 



Rt{x, u) = Wt(x, u) 



(7) 



When a more complex wavefront hits the surface, we have to de- 
compose the incoming wavefront Ri into plane waves and trans- 
form each to reconstruct the outgoing wavefront. This transforma- 
tion is similar to the rendering equation: 



Ro{x,Uo) ^ J Wt{x,Uo,Ui)Ri{x,Ui)dui 



(8) 



where Ro is the transmitted wavefront for an incoming wavefront 
Ri. Wt{x,Uo,Ui) is solely dependent on the micro structures of 



geometry t(x). Under the assumption that the structure is suffi- 
ciently thin, this transformation becomes angle shift-invariant; the 
transmitted pattern simply rotates with the incident rays. Note that 
this property does not necessarily imply isotropy, as the transmitted 
wavefront does not have to be symmetric. An angle-shift invariance 
assumption is valid when the lateral size of a patch is significantly 
larger than the thickness of the microstructure. Under this assump- 
tion the transmission function is given by 



Ro{x, Uo) 



j Wt{x,Uo - Ui)Ri{x,Ui)dui 



(9) 



This transformation is also calculated using the WDF as in Eq [2] 
Note that Wt{x,u) may contain positive as well as negative 
real- valued coefficients. However, after the integration over a 
finite neighborhood, the intensity value always becomes non- 
negative [Bastiaans 1997 1. For reflections we assume that the re- 
flected wavefront is a mirror of the transmitted wavefront, and 
therefore: 



Ro (^X , Uo 



: J Wt{x, -Uo ■ 



- Ui)Ri{x, Ui)dui 



(10) 



This leads to the light transport defining wave rendering equation 
for thin microstructures (compare with Eq. [TJ, 



Ro{x,Uo) ^ Re(x,Uo)-\- J Wt{x, -Uo-Ui)Ri{x,Ui)dUi 



2.1 Converting Microstructures to WBSDF 



(11) 



In order to create the WBSDF, we need to know the exact mi- 
crostructure t{x) of the material. The microstructure is defined by 
an amplitude a{x) and phase $(x) with 



t(x) = a(x)e^^("\ 



(12) 



Equation [12] allows us to calculate the WDF for the microstructure 
by Eq.[2] Finally, the BSDF is calculated using a basis transforma- 
tion as 

p{x, Oi, Oo, X) (x, , (13) 



with positive sin(^o) in the case of transmission and a negative for 
reflection. 



The amplitude a(x) of the microstructure defines transmissivity, 
reflectivity and absorption. For example, an open aperture is 
a function that is 1 within the opening and elsewhere. The 
component $(x) indicates the phase delay introduced due the 
refraction index or the thickness of the material at position x. Fig- 
ure|3ja) illustrates the phase function for a sinusoidal grating in 2D. 

As an example, consider the sinusoidal grating in Figure [3] We can 
formulate the microsurface as 



t(x) — e' 



i^sin(27rf ) 



i27rq - 



(14) 



where m is the maximum height of the grating, p is the period and 
Jq is the Bessel function [Goodman 2005 1. Using Eq. [2]we can 
calculate its Wigner Distribution Function as 

oo oo 

xe--.-^ ^^^s(u-^^y (15) 
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plane 
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Figure 4: A schematic overview of the approximated surface field, 
where ki is the center wave vector incident at the elevation angle 
Oi. h(x) represents the height difference from the target plane to 
the actual surface. 



As a final step we convert this to the Wave BSDF using Eq.fTS] 



2.2 Creating a Statistical WBSDF 

We can easily express the microstructure for a simple surface. 
Eq. (|2| implies that the BSDF can be computed provided the exact 
surface profile of the material. However, the exact microstructure 
may not be known for challenging objects. One way to circum- 
vent this is to represent such surfaces statistically via autocorrela- 
tion and standard deviation (articulating periodicity and roughness, 
respectively). In this situation, we can compute the WBSDF with a 
statistical average as 



/ 



t\x + 



-i27vx' u 



dx\ (16) 



where ( ) denotes average. Depending on the surface properties and 
rendering environments, different types of statistical averaging can 
be used. In general, Gaussian statistics is assumed and statistical 
parameters such as standard deviation a or autocorrelation length 
can be tuned. Typically a standard deviation of height an is nor- 
malized by the wavelength A of light as 



A ' 



(17) 



This statistical average approach has been used in the diffraction 
shader |Stam 1999] and BRDF estimators [ Hoover and Gami7| 
2006|. Note that (t (x + x/2) t* (x - x ' 12)) is mutual intensity 
J{x,x). According to Goodman 11984) , who calculated the sta- 
tistical properties of laser speckles basea on the tangent-plane ap- 
proximation, the surface field can be expressed as 



,/ \ / \ ik^x (l-\-cosOi)h(x) 

t{x) — a(x)e * e ^ ^ w w 
— a(x)e 



ikix ia{x) 



(18) 
(19) 



where ki is the center wave vector incident at the elevation angle 
Oi and h{x) is the height difference between the tangent-plane and 
the actual surface as illustrated in Figure |4] We assume no shadow- 
ing and no multiple interreflections. Also, the roughness should be 
cr > 1. 

From Eq. (T9\ we can derive the phase variance 



[27V a {1 + cos ^2 



and the phase autocorrelation 

\27V_ 
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Ra{x) = 



(1 + COS Oi 



Rh{x) 



(20) 



(21) 



where Rh{x') is the autocorrelation of the surface height (Figure 
[6fa)). If we assume a Gaussian distribution of surface height we 
can now formulate the correlation function as 



7(x ) 



ikix' ^^ict(x) — ct(x — x') 



where 



_ ^ikix' ^-a^[l-ph(x') 



, Rh{x') 

Ph{x ) = ^2 • 



(22) 
(23) 

(24) 



Since the WDF is defined with respect to the correlation function 
as in Eq. (|2]), we can derive the WBSDF from the WDF of the cor- 
relation function as 



-i27ru-x' 



exp 



(25) 



A 



2.3 Benefits 



Near and far field: The WDF as w ell as the WBSD F are 
conserved along rays in the paraxial region |Bastiaans 2009| , hence 
they are valid in both the near-field and far-field. In the far-field, the 
observed wave at a single point is only dependent on angle, and is 
essentially independent of the distance from the grating. Near-field 
corresponds to the Fresnel region in optics (not to be confused with 
the near-zone in optics where the evanescent field is still strong). 
The near-field is the region close to the grating, where the wave's 
distance to the grating also influences the observed pattern. In 
contrast to previous diffractive BRDFs [Stam 199 9], the WBSDF 
does not require an assumption that the object and receiver are at in- 
finity. Near field effects, such as lenses, require more computation 
than far field effects. Hence, importance sampling is highly desired. 




\ 



Figure 5: Renderings using the WBSDF under global illumination. 
(Left) Light reflected off a CD on a diffuse wall, (Right) light emitted 
from an area light ( top) is reflected off the left CD onto the right CD 
and then to the floor. 



polarized light. The current model of the WDF only encodes spa- 
tially varying phase and does not take temporal phase into account. 
To model transient responses, one can incorporate the spatial phase 
from the WDF with temporal frequency. 

Using the WBSDF, a single reflected ray off a surface may con- 
tain negative radiance, which is physically impossible. However, 
by integrating over a finite neighbourhood, for example in a camera 
aperture, the total amount of received light always becomes non- 
negative due to the properties of the WDF [Bastiaans 2009]. There- 
fore, our model does not support pinhole camera models for direct 
lighting. More details on this topic are presented in the supplemen- 
tal material. 

3 Wave-Based BRDF in Practice 

We show representative ray-based wave effect rendering as exam- 
ples for global illumination and sound propagation. 



Importance sampling and global illumination: The WBSDF 
uses instantaneous diffraction; i.e., the magnitude of energy 
scattered in all directions is known at the surface. This allows 
importance sampling for an efficient light simulation and is also 
suitable for rendering global illumination. 

Light coherence: Natural scenes are composed of incoherent 
sources. This means that each single point source (at a reason- 
able distance) is coherent with itself and mutually incoherent with 
other point sources. Simulating interference effects normally re- 
quires rendering the scene for each light source independently, and 
then summing up images to create the final result. Because creating 
diffractive effects for a single light source only consists of summing 
up photons (positive or negative), our technique can simultaneously 
render coherent light as well as multiple incoherent light sources, 
such as area lights. In rare scenarios, multiple sources are mutually 
coherent. Here, we treat all the sources jointly to compute a single 
emissive WBSDF 

2.4 Limitations 

We use the paraxial approximation, where incoming and scattered 
(diffracted) light propagate not far from the optical axis. To im- 
prove the accuracy in the non-paraxial zone, one can use the angle- 
impact WDF I Alonso 2004 1. Regarding polarization, we only con- 
sider linearly polarized light. Adopting the coherency matrix | Tan- 
[nenbaum et al. 1994J , we can extend our method beyond linearly 



3.1 Light Rendering 

Diffraction and interference: The WBSDF indirectly retains the 
phase information after diffraction, and defers the interference com- 
putation. This strategy supports interference after global illumina- 
tion naturally. Figure [T] shows the diffracted wavefront from a CD 
undergoing refraction. As we describe later, instantaneous diffrac- 
tion enables importance sampling. 

Figure[5]shows 2 CDs interacting with each other and illuminating a 
wall. We model the CD as a phase grating with a pitch of 2 /xm. We 
believe this is the first practical demonstration of global illumina- 
tion of wave phenomenon, which is achieved using an unmodified 
PBRT. 

Figure [6] demostrates the statistical approach of the BSDF. In the 
autocorrelation of the surface we used the function 

Rh = -(x/a)^ (26) 

as shown in Figure |6ja). The BSDF is calculated for Oin and Oout 
(Figure [6fb)) and mapped onto spheres, each with a different stan- 
dard deviation. 

Rainbow holograms: As we can create the WBSDF of any 
diffractive optical element based on its microstructure, we can de- 
rive a WBSDF for a holographic surface as well. As diffraction 
is computed instantaneously for the entire surface, we do not keep 
track of travel distance from a given point to the hologram surface, 





Figure 6: Rendered results using the statistical approach of the 
WBSDF. (a) Example Ru function describing the basic surface 
height of elements to reconstruct the micro structure, (b) The sta- 
tistical WBSDF of element (a) plotted in Oi and Oo for a specific 
standard deviation, (c) Rendered result of spheres on which this 
BRDF is mapped. Spheres 1 to 4 have respectively a standard de- 
viation of 10, 8, 6 and 4. 



as OPD does. The hologram can be rendered directly from the WB- 
SDF and creates interference inside the aperture of the camera lens. 
Note that t here are many different w ays of computing the WDF of 
the object [Plesniak and Halle 2005 1. For simplicity, here we en- 
coded incoherent objects in the hologram; however any arbitrary 
wavefront information can be stored in the WBSDF. 



and point light positions. We then applied the PSFs to a computer 
generated image by convolution, simulating the appearance of the 
scene when viewed with this lens. Note that in this example, no 
chromatic aberration was assumed. Hence, the color dispersion 
results solely from diffraction. Although, we only consider thin 
lenses in this particular example, we can easily simulate a series of 
thick lenses as well. The global illumination ray-tracing takes care 
of refraction at air-glass interfaces and diffraction by an aperture is 
included by the WBSDF. 
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Figure 8: Near field effects in camera lenses show the effect of an 
aperture creating diffraction and chromatic dispersion. (Left) PSFs 
of F/5.6 lens for lights at two different depths. (Right) PSFs ofF/11 
lens for two different depths. 

4 Implementation 



3.2 Simulating Camera Optics 

In camera optics, aberrations as well as diffraction affect the Point 
Spread Function (PSF) of the system. Hence, conventional optics 
design software suites provide analytical tools for aberrations and 
diffraction. However, they treat the two separately; i.e., ray tracing 
uses spot diagrams as an estimate of the PSF, whereas diffraction 
models produce PSFs directly from Fourier optics analysis. Note 
that a Fourier optics analysis typically assumes a shift invariant PSF, 
hence the diffraction for an off-axis PSF is assumed to be identical 
to the on-axis PSF, which is not true in practice. 

In addition to Oh et. al. |2010|, our approach can simulate both 
aberrations and diffraction simultaneously for any location of the 
image plane, since it is based on ray-representation and able to 
include diffraction. We derived the WBSDF from the geometric 
structure of a camera lens with a circular aperture and calculated 
the PSF. Figure [8] shows spatially varying PSFs for different lenses 




Figure 7: Rainbow hologram renderings from different viewpoints. 
Instantaneous diffraction using the WBSDF pre-computes the view- 
dependent appearance. 



We achieve global illumination effects using an unmodified PBRT 
framework [Pharr and Humphreys 2004 1. The only minor change 
is that we comment out the check for non-negative radiance. To 
generate specific diffractive material, we simply created a new ma- 
terial plugin for PBRT. 

We used photon mapping for Monte Carlo simulation of global il- 
lumination. Due to the nature of this technique, photon mapping 
performs well with our BSDF. Photons can convey negative energy 
to achieve destructive interference. WDF formulation ensures that 
non-negative final energy is gathered at a single point. This non- 
negativeness also applies to the gathering of photons in the photon- 
map. 

For materials with a separable WBSDF for x and y coordinates, we 
precomputed the WBSDF and saved it as a 2D look-up table, which 
allows fast calculation, importance sampling (Figure [9|, and real- 
time rendering of single bounce diffraction on the GPU (Figure[TO|. 
As diffration has to be calculated on arrival, we need to sample over 
different locations for each pixel and integrate the values. Mipmap- 
ping calculates these interference e ffects efficien tly. Similar to the 
diffraction shader implementation |Stam 1999| , the colorf ull 0th 
order highlight is added with a simple anisotropic shader f Ward| 
|1992| . Table |2] gives an overview of rendertime used for generat- 
ing the images and audio examples. The pdf is calculated using the 
absolute values of the intensity in the BRDF lookup table. 

Eq. |T3] ) provides a function in position, angle, and wavelength. 
For structured surfaces, like sinusoidal gratings, we have a closed 
form solution to find which wavelengths have non-zero intensities. 
For computational efficiency, we tabulated the WBSDFs in terms 
of position and angle, and reduced the wavelengths to RGB values 
according to the camera response curves. The lookup tables, used 
for creating the examples in this paper, have a sampling rate of 
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Renderingtme 


5100 s 
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Table 2: Rendering time and number of photons used for each 
scene. Implementation is done using PERT. All these scenes are 
rendered single threaded on a 3 GHz core. The GPU example was 
done on a NVidia 8400M with 640x480 resolution. 
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(a) Without importance sampling (b) With importance sampHng 

Figure 9: Benefits of importance sampling (a) Without importance 
sampling, even after 14,000 s result does not converge (b) With, 
1,300 s. 






Figure 10: CD rendered using the WBSDF in real-time on a GPU 
using look-up tables and mipmapping. 



0.02° in angular domain and 48 nm in the spatial domain. For 
non-structured surfaces such as rainbow holograms, we computed 
Eq. ( pj] ) over 30 wavelengths and stored only RGB values using 
the appropriate weights. 



Validation We compared two representative diffraction examples 
with the one computed by Fourier optics and verified the accuracy 
of our technique. We used photon mapping to simulate Young's 
double slit experiment and light diffraction from a rectangular aper- 
ture (Fraunhofer diffraction). The absolute intensity error, mainly 
due to limited numerical precision, between the Fraunhofer diffrac- 
tion pattern and the rendering with the WBSDF is less than 1%. 
These results are presented in Figure pT] 

5 Comparison 

5.1 Comparison with OPD 

Even though optical path difference (OPD) provides accurate re- 
sults of interference, it requires significant modifications to tra- 
ditional rendering systems and is not able to perform importance 
sampling. Traditional raytracers commonly do not take distance 
or phase into account, and in order to implement OPD we need to 
make changes to the framework. The precision required for path 



length (for each wavelength) becomes challenging. In addition, 
OPD cannot exploit importance sampling. Determining the propa- 
gation direction of rays with dominant intensity in the presence of 
diffraction is challenging. OPD-based techniques must uniformly 
sample all the outgoing directions. Consider the example of the si- 
nusoidal gratings of Figure [3] The WBSDF reveals the few impor- 
tant directions using our instantaneous diffraction theory. In more 
complex, global illumination settings, WBSDF-based approach is 
faster compared to OPD by the same factor as shown in Figure |9] 

5.2 Comparison with Diffraction Shaders 

Diffraction shaders (DS) pioneered the idea of fast and practica l 
rendering of the far field diffraction of a single bounce |St am 1999) . 
It turns out that the DS approach is a special case of our method 
where the light source and the observer are at infinity. DS converts 
parallel ray bundles predicted by the WBSDF and treats them as a 
single ray beam (Figure [T2|. The WBSDF maintains a higher res- 
olution representation and hence we can apply different integration 
kernels, transforms for propagation, and scattering models. The 
WBSDF also provides flexibility depending on incident/outgoing 
rays, camera geometry, computational speed and tolerance of error. 
The WBSDF supports a statistical model similar to the one used 
in DS for phase variations. Based on the assumption that the light 
source a nd the obse rver are at infinity, we use the same notation as 
in Ref. | ]Stam 1999| and present the relationship with the WBSDF 
as follows: 

/ [ku) = ip {ku) V^* {ku) 

^11' 



^ikwh(p-\- ^) ^-ikwh(p- ^) ^iku q ^^^^ 



I 



Wd{x,u)dx, 



(27) 



where / is intensity, ip is field, k = 27v/X, u' = sin^i — sin ^2, 
w = — cos^i — cos^2,P — {x + x')/2, q — x — x', and Wcz(x, u') 
is the WDF of q^^^^^"") with respect to sin Oi - sin O2, which rep- 
resents the reflectance of the surface. In our formulation, as men- 
tioned in Sec. 2, the outgoing WDF is written as 



Ro (x, Uo) 



j Wt {x,Uo]Ui) 



R{x,Ui)dui. (28) 



Assuming that a plane wave is incident on a surface (Ri{x, Ui) — 
5 (ui — sin Oi /A)) and the observer is at infinity as in the diffraction 
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Figure 11: Validation of the WBSDF (a) Visual comparison of 
photo and rendering of diffraction due to a laser through a rect- 
angular aperture, (b) Rendering of (a) and its ID plot, (c) Two 
slit experiment and ID plot. The ID plots are compared with the 
Fraunhofer diffraction and have an absolute intensity error of less 
than 1%. 
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Figure 12: Diffraction shaders is a special case of the WBSDF. (a- 
b) It pre-encodes the result of far-field interference into energy of 
the outgoing ray. These rays cannot destructively interfere anymore 
in a later stage. ( c-d) WBSDF indirectly encodes phase information 
into the outgoing ray by introducing potentially negative radiance, 
allowing the rays to interfere later for global illumination. 



shader equations {Uo — sin Oo/X), we obtain the reflected light as 



I(uo)^ j Ro{x,Uo)dx ^ J 



^^1^ sm^.smt^i^^ (29) 



In the implementation of diffraction from CDs, we assumed angle- 
shift invariance, where the phase function ip(x) does not depend 
on the incident angle. If we use the same setup assumptions as DS 
(i.e., a plane wave incident on a camera at sufficiently far distance), 
then our result is comparable to the one generated by DS, especially 
in the paraxial region. 

The WBSDF analysis also points to new ways to improve the DS 
approach. DS can now be used for reflection using Eqs. |9]) and 
|To] ). In the examples generated by DS, height-maps were used to 
model phase delays for pre-computing interference from parallel 
light rays, ignoring the amplitude component of the micro structure. 
The amplitude controls the amount of light that is absorbed or scat- 
tered. To account for the amplitude variation, we can add an a(x) 
term in Eq. p2] ). Despite these extensions, DS remains applicable 
only in the far field, since the integration is already performed. DS 
is also less suitable for simulating audio diffraction (large wave- 
length) or simulating a PSF of a camera's optics (larger integrating 
cone for receiving patch). 

5.3 Comparison with Augmented Light Fields 



Augmented light fields |Oh et al. 20T0) (ALF) is the first model 
to use the WDF for rendering. It uses a simplistic two plane light 
field parametrization: one for diffraction grating plane and one for 
receiver plane. Hence, it is limited to demonstrations of a planar 
wavefront transmitted via the first plane of a diffractive occluder. 
They use a 'destination based' approach (i.e. differed diffraction) 
using an OpenGL fragment shader. The ALF is computed in a back- 
ward manner: at each point on the receiver, the shader computes the 
incident ray-bundle (slice of the lightfield) using an analytical for- 
mula for the diffraction pattern from the source towards the receiver. 
This differed diffraction is very similar to OPD, so it provides no 
importance sampling and the paths to trace grow exponentially with 
each bounce. The ALF work suggests photon mapping as a possi- 
ble extension, but it is treated as a procedure that would modify the 
scatter stage of a renderer. By encoding a new BSDF, we make no 
change to the renderer. Instead of a new rendering strategy, WB- 
SDF encodes the micro structure into a new reflectance function in- 
dependent of other elements or illumination in the scene. This al- 
lows a wide range of effects not shown with ALF such as reflection, 



emission, multi-bounce, importance sampling, sound rendering and 
a natural fit with the rendering equation. Refer to Table 1 for more 
details. 

5.4 Comparison with Edge-diffraction 

Edge diffraction is a rudimentary form of importance sampling used 
in audio rendering. It can be improved by using the WBSDF as 
shown in Figure [13] This becomes even more significant when 
phase is involved. Edge diffraction techniques underestimate mid- 
air diffraction and allow strong ray bending at edges, whereas our 
technique creates diffraction through entire features, ensuring accu- 
rate diffraction results. Edge diffraction cannot achieve the unusual 
effect of the sound of a door closing. Our current method does not 
take traveling distance into account and as a result does not exhibit 
damping. 
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Figure 13: WBSDF is more accurate than edge diffraction, (a) 
Edge diffraction only diffracts rays at corners and (b) ignores mid- 
air diffraction, (c-d) WBSDF creates midair diffraction. 



6 Conclusion 

We describe a new representation of BSDF that greatly simplifies 
simulations of wave-phenomena in ray-based renderers. We bring 
the wave phenomena into the realm of the rendering-equation which 
supports global illumination, including benefits such as importance 
sampling. Additionally, we provide a detailed comparison of many 
wave phenomena renderers. We feel that we have methodically in- 
vestigated the rendering of wave phenomena by proposing an easy 
and efficient solution, demonstrated examples in multiple domains 
(light/sound), investigated and solved implementation issues, and 
mathematically expressed out work's relationship with diffraction 
shaders. 

We see many promising directions future exploration. The WBSDF 
can also be used for other areas of spectrum such as microwaves 
and x-rays or other wave-related disciplines, including the tran- 
sient response for fluid wave diffraction. The simulation of camera 
optics can be extended to support a full camera model in existing 
rendering frameworks. Sound rendering extensions include atten- 
uation, reverberation, non-planar area sources, Doppler effect and 
modeling of head-related transfer function (HRTF) using our for- 
mulation for easy ray-based analysis. Advanced surface scattering 
models, e.g., tangent-plane approximation or phase-screen approx- 
imation, will potentially extend the applications of the WBSDF and 
improve the accuracy to the physical models. Although current re- 
flectance field scanning methods use a ray-based understanding, we 
hope WBSDF methods will inspire novel capture and inverse prob- 
lems in analysis of real world objects exhibiting wave phenomena. 
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Appendix 

Wigner distribution function for ID element t{x),in Matlab. 

function W = WDF(g) 
N = length(t); 

X = (((0:N-l)-N/2)*2*pi/(N-l)); % Generate linear vector for shift 
X= (0:N-l)'-N/2; 

Gl = ifft( (fft(g)*ones(N,l)).*exp( i*x*X/2 )); % create g( x + x/2' ) 
G2 = ifft( (fft(g)*ones(N,l)).*exp( -i*x*X/2 )); % create g( x - x/2' ) 

W = fft(Gl.*conj(G2), [], 1); % calculate WDF 



